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We calculate the conductances of a three-terminal junction set-up of spinless Luttinger liquid 
wires threaded by a magnetic flux, allowing for different interaction strength g^ ^ g in the third 
wire. We employ the fermionic representation in the scattering state picture, allowing for a direct 
calculation of the linear response conductances, without the need of introducing contact resistances 
at the connection points to the outer ideal leads. The matrix of conductances is parametrized by 
three variables. For these we derive coupled renormalization group (RG) equations, by summing up 
infinite classes of contributions in perturbation theory. The resulting general structure of the RG 
equations may be employed to describe junctions with an arbitrary number of wires and arbitrary 
interaction strength in each wire. The fixed point structure of these equations (for the chiral Y- 
junction) is analyzed in detail. For repulsive interaction [g, gs > 0) there is only one stable fixed 
point, corresponding to the complete separation of the wires. For attractive interaction {g < 
and/or (73 < 0) four fixed points are found, the stability of which depends on the interaction 
strength. We confirm our previous weak-coupling result of lines of fixed points for special values of 
the interaction parameters reaching into the strong coupling domain. We find new fixed points not 
discussed before, even at the symmetric line <? = gs, at variance with the results of Oshikawa et al. 
The pair tunneling phenomenon conjectured by the latter authors is not found by us. 
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I. INTRODUCTION 

Electric circuits of one-dimensional quantum wires are 
described by laws quite different from those of corre- 
sponding classical systems. These laws are non-local, of- 
ten quantized, and are generally difficult to derive. An 
important first step is the exploration of the properties 
of junctions of wires. The linear response conductance 
of a (symmetric) two-terminal junction in the limit of 
zero temperature has been found to be either zero (re- 
pulsive) or unity (att ract ive interaction) in units of the 
conductance quanturrp321. Most of the existing theories 
employ the exactly solvable Tomonaga-Luttinger liquid 
(TLL) model, in which backward scattering is neglectecP. 
Moreover, for simplicity the fermions are assumed to be 
spinless, in which case backward scattering only gives 
rise to a renormalization of the forward scattering am- 
plitude. Furthermore, most of the early works on the 
problem of two TLL- wires connected by a junction em- 
ployed the method of bosonization, or, for special values 
of the interaction pa ram eter, used a mapping on to ex- 
actly solvable models^. 

An alternative formulation using a purely fermionic 
representation in terms of scattering states and a renor- 
malization group (RG) approach for the conductance has 
been pioneered by^j. The latter formulation has several 
advantages over the bosonic approach. Perhaps the most 
important difference is that it allows to describe the phys- 
ically relevant system of a quantum wire (with interac- 
tion between the fermions) adiabatically coupled to reser- 
voirs (with negligible interaction), whereas the bosonic 
approach is formulated for the infinite system. The effect 



of a finite wire length has only been convincingly derived 
for a clean wire ^^ ; in the case of wires connected by 
a junction a plausible, but unproven ad hoc procedure 
has been proposed^*' . As will be discussed below we have 
reasons to suspect that the latter procedure is not always 
correct. The fermionic scattering state approach is gen- 
erally applicable, allowing for the st udy o f wires with two 
impuritiesIH^ of multi-lead junctionJi^El! as well as the in- 
clusion of spin, of backward scattering,'^ of out of equilib- 
rium situationSjli^ point contact in a quantum spin Hall 
insulatopi^. The original model calculatiorP has been 
limited to weak coupling, ft may be shown, however, 
that the method can be extended to work in the strong 
coupling regime as well, by summing up an infinite class 
of leading contributions in perturbation theory^®. In all 
cases considered so far it has been found that the method 
provides results in agreement with any available exact re- 
sults, which builds confidence in its validity. As shown 
iiF^, in the case of a two-lead junction, the contributions 
neglected in the above-mentioned partial summation are 
subleading and may be shown to vanish when the fixed 
points are approached. 

Generally, the transport behavior of one-dimensional 
systems at low energy/temperature is dominated by only 
a few fixed points (FPs) of the RG flow. The most com- 
mon and intuitively plausible FPs are those with quan- 
tized conduction values, G = 0,1, in units of the con- 
ductance quantum Gq — e'^/h (we consider systems of 
spin- less fermions). There may, however, appear addi- 
tional FPs associated with noninteger conductance val- 
ues, usually in the case of attractive interaction. A par- 
ticularly well-studied case is that of a symmetric three- 



lead junction with broken chiral symmetry, as induced 
by a magnetic flux threading the junction. This latter 
problem has been studied by the method of bosonization 
i] jio | i7 |^ and by the functional RG method in^^. These 
authors have identified a number of new FPs, for attrac- 
tive interaction. The discovery of these new FPs has 
sparked interest in the problem of mapping out the com- 
plete fixed point structure of the theory, even though 
most of the interesting new physics appears to be in the 
physically less accessible domain of strongly attractive 
interactions. Only the fully symmetric situation of iden- 
tical wires {g = gs, see below) has been considered irPSESl. 
As we have shown iniSI for weak coupling and will show 
below for any coupling strength, important new physics is 
missed by restricting the consideration to the symmetric 
case. 

One result found by"'^ deserves special attention: it 
is the claim that a new fixed point named D (Dirich- 
let) exists for strongly attractive interaction (Luttinger 
parameter X > 3), which is interpreted as describing 
multiparticle transport. As a consequence, the conduc- 
tance values at this FP are outside the domain of values 
allowed by a single particle S-matrix description. Put dif- 
ferently, this requires the excitation of particle-hole pairs 
at the junction, with say the particle moving into lead 2, 
while the hole is going into lead 3. In our view, at least at 
zero temperature, and in the linear response regime, such 
processes are not occurring, because the phase space for 
particle-hole excitations tends to zero in this limit. As 
we will show below, we do not find a fixed point with 
the characteristics of D, but in contrast we find two new 
fixed points, not discussed in Refs. [10120) . 

In principle, the Y-junction set-up may be realized ex- 
perimentally by a one-dimensional tunneling tip contact- 
ing a quantum wire. A magnetic flux threading the junc- 
tion may be created by local magnetic moments at the 
junction, although this has not been realized experimen- 
tally so far. 

In this paper we report results of an RG treatment of 
electron transport in the linear response regime through 
a junction of three spinless TLL wires threaded by a mag- 
netic flux. We employ a fermionic representation as de- 
scribed indetail irfI3. The method has formerly been 
used irflZMl [j^ ^j^g ^g^gg Qf Y-junctions. More recently 
this approach has been used to derive the RG equations 
for the conductances of a Y-junction connecting three 
TLL wires in the absence of magnetic flux, for weak 
interaction^, and in a recent work for any strength of 
interaction^^. In^^ it was found that even in weak cou- 
pling, but beyond lowest order, interesting new struc- 
tures appear. Even a weak higher order contribution may 
change the RG-flow in a dramatic way. In the present 
paper we describe a similar effect: an asymmetry of the 
three wires of a chiral Y-junction to the effect that in the 
tip wire the interaction strength g^ is different from that 
in the main wire, g, allows to access certain regions in 
interaction parameter space where a whole line of fixed 
points rather than a single fixed point is stable. This hap- 



pens for attractive interaction only {g, g^ < 0) The line 
of stable fixed points is connecting two fixed points at two 
special manifolds of interaction values, (a) g = 0; (?3 < 0, 
or (b) whenever the condition K^^ = 2 — {K + ii'~^)/2, 
1 < K < 2 \s met. Along these fixed point lines the 
FP values of the conductances are continuously vary- 
ing. Results on this new aspect of transport through Y- 
junctions within weak coupling have been reported iiPSl. 
In the present work we extend our theory into the strong 
coupling regime. We derive a general expression for the 
scale-dependent contribution to the conductances, valid 
for a junction with any number of leads and any values 
of interaction gj in lead j . This result allows to derive 
RG-equations for the set of independent conductances 
in the given case. We go on to evaluate these general 
expressions for the three-lead junction of the K-type (in- 
teraction g in the main wire and 93 in the tip wire) in 
the presence of magnetic flux. Then we determine the 
six physically allowed FPs, and discuss their stability as 
a function of position in coupling constant space. 

In order to compare our results with the recent flnd- 
ings oPSl for a fully asymmetric junction, we extend our 
analysis by allowing for full asymmetry of the S-matrix 
at the junction, while keeping equal interaction in the 
main wires. This minimal modification enables us to elu- 
cidate the differences between our conclusions and those 
iipSl. We show that the set of RG equations for the most 
general asymmetric chiral case can be easily derived in 
our formalism, and allows us to discuss the robustness 
of our results obtained for the partially asymmetric case 
(51 = 52) with respect to further asymmetry. A full anal- 
ysis of the case gi ^ g2 is, however, beyond the scope of 
the present study. 



II. PERTURBATION THEORY FOR 
CONDUCTANCES 

A. The model 

We consider a model of interacting spinless fermions 
describing three quantum wires connected at a single 
junction by tunneling amplitudes in the presence of a 
magnetic fiux piercing the junction. In the continuum 
limit, linearizing the spectrum at the Fermi energy and 
including forward scattering interaction of strength gj 
in wire j, we may write the TLL Hamiltonian in the 
representation of incoming and outgoing waves in lead j 



(fermion operators '4>j.im '4'j,out) ^s 
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i/™* = 27r{g[pji^ + p2%] + g3P3P3}<S>{x; 



(1) 



~L,-l). 



Here *»„ = (V'l.m, V'2,m,^3,m) denotes a vector oper- 
ator of incoming fermions and the corresponding vector 



of outgoing fermions is expressed through the S'-matrix 
as ^^outix) = S ■ \l/i„(— a;) . In the chiral representation 
we are using positions on the negative (positive) semi- 
axis corresponding to incoming (outgoing) waves. We 
consider quantum wires of finite length L, contacted by 
reservoirs. The transition from wire to reservoir is as- 
sumed to be adiabatic (i.e. produces no additional poten- 
tial scattering). The junction is assumed to have micro- 
scopic extension I of the order of the Fermi wave length. 
Inside the junction interaction effects are neglected. This 
is expressed by the window function Q(x; —L, —I) = 1, 
if — L < X < —I, and zero otherwise. The regions 
X < —L arc thus regarded as reservoirs or leads la- 
beled j = 1,2,3. We denote the interaction constants 
9i = 92 = 9 from now on, and put the Fermi velocity 
vp = 1. The various incoming and outgoing channels are 
illustrated in Fig. 1 elsewhere.'^^'^SI The interaction term 
of the Hamiltonian is expressed in terms of density op- 
erators pj,in = *+Pj* = pj, and pj^out = VP+Pj* = pj, 
where pj = S^ ' Pj ' ^ and the density matrices are given 
by {Pj)ai3 = Sa^Saj and (p,)q/3 = S'^jSjp. The S'-matrix 
characterizes the scattering at the junction and (up to 
irrelevant phase factors) has the structure (se^^-) 



S = \ t21 n t23 

^31 ^32 r2 



(2) 



A convenient representation of 3x3-matrices is in 
terms of Gell-Mann matrices Xj, j — 0,1,..., 8, the 
generators of SU{3) (see22l). Notice that the interac- 
tion operator only involves A3 and As (besides the unit 
operator (Ao)q^ — y^2/3Sap)- We note Tr(Aj) = 0, 
Tr(Aj Afc) = 26jk, j = 1, . . . , 8, and [A3, As] = 0. We intro- 
duce a compact notation by defining a three-component 
vector A = (A3,As,Ao), in terms of which the densities 
may be expressed as pj = 
3x3 matrix R is defined as 



l/^Su-^iM-^A" where the 
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(3) 



and has the properties R ^ = R"^, deiR — 1. The 
outgoing amplitudes are expressed in the analogous form 



with A, replaced by A,- 



S+ ■ Aj • S. 



With the aid of the 



\j the S'-matrix may be parametrized by eight angular 
variables (see^S and Appendix [A| . For the case under 
consideration only three of these, 6,tlj,£,, are relevant: 
S = gjA2j/2ga3(77-v)g»A5egjA2(77-0/2^ Tj^g corresponding 

elements of the S-matrix are given by 



ri = - (cos 6' + e"'' ) sin ^ , r2 — cos 6 , 
*13 = ^32 = » cos - sin 6* , ^23 = ^31 = «sin-sin6l, (4) 
ti2 = cos e cos^ - - e"'' sin^ - , ^21= ^12 le^^-^ ■ 



As will be shown below, the S-matrix and therefore the 
angular variables 9,ip,^, will be renormalized by the in- 
teraction. 



B. Linear response conductances 

In the linear response regime, the conductances 



Cjk 



dy {(Pji-L) - pJL))pk{y))u^^Q 



relate the current /-,- in lead j (flowing towards the junc- 
tion) to the electrical potential T4 in lead k, Ij = CjtVk- 
Here the combination (j)j{~L) — Pj{L)) means that the 
current is measured at x = —L, and the integration over 
y for the incoming densities is over the length of the ideal 
leads, to which the electric potential T4 is applied. De- 
noting the frequency of the applied external field by uj, 
we consider the d.c. limit toL — >■ 0, when the value of 
L does not play a role-t^^ The quantum averaging (. . .) 
involves taking the trace over the matrix products. The 
conductance matrix C has only three independent ele- 
ments, Gaa = 5(1-0), Gbb = §(1-^) , and Gab = c/^S, 
relating the currents /a = ^(/i — 12), lb — |(/i+/2 — 2/3) 
to the bias voltages T4 = (V1-V2), Vb = 5(^1-^^2-2^3). 
In compact notation we may define a 3 x 3 matrix of con- 
ductances 



G = 



Gaa Gab 

— Gab Gbb 





i(l 



V3 






(5) 
The connection between matrices C and G comes from 
the observation, that the nonzero elements of the matrix 
C^ = R-^CR are essentially the reduced conductances: 



Cji — 2Ga 



^22 



\Gbb , C«2 = -^2^1 = A/3Gab , where 



the numerical factors arise due to the physically moti- 
vated asymmetric definitions of the currents and voltages. 
It may be shown that the parametrization of the conduc- 
tance tensor in terms of a, 6, c follows quite naturally, by 
observing that the initial conductances following from 
the Kubo formula are given by Gju = Sjk — Tr(p ''p^) = 
^jk — \S,fk\'^ where the superscript r denotes that the quan- 
tity is fully renormalized by interactions^S. By express- 
ing Yjk = Tr(p/pfc) in terms of Y^^ = iTr(A;;A,) as 



Y = (R • Y • R"^), we find by comparison with the con- 
ductance matrix, C'^ = 1 — Y , that Y-'^ has nonzero el- 
ements given by the conductance parameters introduced 
above: 



Y^ 




(6) 



From now on we wiU drop the superscript r , with the un- 
derstanding that all quantities are renormalized. There- 
fore, we may use Y^ to represent the conductances in 
the renormalization group analysis below. 

By substituting the iS-matrix in the form (pi), Q into 
the definition of Y^, we find the relation of the a, &, c 
to the Euler angle variables: a — — ^(cos^ 9 + 1) cos-^ ^ + 

1),C: 



,V3, 



I COS 
2 



0COS^. 



COS0COS V'sin^ ^, 6 = i(3cos^ >^ ^j,^— ^ 
We find therefore that a, b, c are confined within the re- 
gion a e [-1,1], b e [-1/2,1], c e [-73/2,^3/2]. The 
physically allowed points in a-6-c-space, satisfying the 
conditions |cos0| < 1, [cosi/'l < 1 and |cos^| < 1, he 
inside a body labelled B as shown in Fig. [ij We may 
trace the above restriction on the allowed values of the 
conductance back to our consideration of energy conserv- 
ing scattering. In the linear response regime and at zero 
temperature only lead electrons at the Fermi level will 
be transported through the interacting wire, and particle- 
hole excitations are excluded. The restriction of the set of 
conductances to the domain of allowed values has impor- 
tant consequences for the stability of certain fixed points 
of the RG flow, as will be discussed below. We mention 
now and will show later that all fixed points are located 
on the surface of the body B in Fig. [l] 

For the isotropic case discussed below we have a = b. 
One possible parametrization of the S-matrix for this case 
is given by Eq. (Bl) The corresponding subspace of al- 



lowed conductance values is the twodimensional domain 
A enclosed by a deltoid curve, see Eq. ( B8 1 and Fig. [2] 



III. RENORMALIZATION GROUP EQUATIONS 

The renormalization of the conductances by the in- 
teraction is determined by first calculating the correc- 
tion terms in each order of perturbation theory. We are 
in particular interested in the scale-dependent contribu- 
tions proportional to A = ln(i/Z), where L and I are two 
lengths, characterizing the interaction region in the wires 
(see above). In lowest order in the interaction the scale 
dependent contribution to the conductances is given byl^l 



1.771 



ffm/A , (7) 



where Cjfc|g=o = S^k-Yjk, the Wjk = [Pj.Pk] are a set of 
nine 3x3 matrices (products of VF's are matrix products). 




FIG. 1: (Color online) Allowed values of conductances, a, b, c, 
lie inside the depicted body, B. The location of the RG fixed 
points A'^, A, X*, M is also shown. The line of fixed points 
connecting FPs A and x'^ is indicated in yellow. 
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FIG. 2; (Color online) Allowed values of conductance a, c 
for the isotropic case. The S- matrices (Bll characterized by 
par ame ters F, </> lead to two quantities a and c, as described 
by (B9l. We took 500 pairs F, (^ at random, these are blue 



points on the plot. As discussed in Appendix pi the domain 
A is bounded by a deltoid curve. 



9ml — gm^mi IS the matrix of interaction constants and 
the trace operation Tr is defined with respect to the 3x3 
matrix space of VF's. As argued irfl^the leading terms in 
each order of perturbation theory are certain diagrams of 
the ladder type, which may be summed up analytically. 
The result is a renormalized interaction matrix g'-'^'^d.^^ 
replacing the bare interaction matrix g in Eq.([7|. The 
components of g'-'^d.der ^^^.g obtained from the following 
matrix equation (see appendix O) 



gladder 



2(Q 



(8) 



We observe that the effective interaction is found to de- 
pend on the conductance components Cjk — 6jk — Yjk 
in a highly nonhnear way. The matrix Q characterizes 
the interaction strength and depends on the Luttinger 
parameters Kj — [(1 — gj)/{l + ffj)]^^^ as 



Qjk = qjS.k, qj = (1 + K,)/{1 - K,) 



(9) 



(Notice that \qj\ > 1 for any interaction strength.) In or- 
der to remove the redundancy in the conductance matrix 
C we now multiply with R^ from the left and R from 
the right to get the components of Y^ in the form 



Yfk = ^f |,.o - ^ E Tr [W.'lw/i] sLr^'-'^A . (10) 
Differentiating these results with respect to A (and 



then putting A = 0) we find the RG equations up to 
infinite order in the interaction 



dA ^^ 



Y^tAw^^wI 



ladder, R 
9ml 



(11) 



Here the WK = {R'^ • W • R}jk are a set of nine 3 x 
3 matrices and the trace operation Tr is defined with 
respect to that matrix space, whereas g^^ '^^' — {R"^ • 
gladder . R^|^^ ^ 2{(Q« - Y^)-^}^i and Q«, Y« are 
scalars with respect to this space. The nine matrices 
WR. are best evaluated with the aid of computer algebra, 
inserting the S'-matrix in terms of the quantities a, 6, c as 
given above (see also Appendix Id|. As a result one finds 
the following set of RG-equations 



da 

dA 

db 

dA 
dc 
dA 



D~'^{a\3b - 1 - 2Qi) + a(3c2 + q{l - b)) + (Qj - 6)(1 + b) + c\2 + q)} 
D-'^{q{l - b){l + 25) + a{b - 1)(1 + 36 - Qi) + c^{2 + 36 + Qi)} , 



-cD-^{l 



}i - 3c^ + 2b{q - 1) + a(2Qi - 36 - 2)} , 



(12) 



where D 



q){b-Qi 



and 



Qi = 



3gg3 -q-'^qa 
2g + 93 - 3 



The RG-equations describe the flow of the conductances 
upon increasing A until a stable fixed point is reached. 
The fixed points are found by putting the right hand 
sides (the /3-functions) simultaneously equal to zero. The 
approach towards any given fixed point is characterized 
by power laws. The fixed point pattern and the power 
law exponents depend on the interaction strength. We 
note that the /3-functions (the right hand sides of Eqs. 
( 12 1 ) are highly nonlinear functions of a, 6, c and of q, q^, 



giving rise to a rich manifold of fixed points and RG flow 
patterns. 



IV. FIXED POINTS OF RG FLOW 

A. Chiral isotropic case 

In the isotropic case, meaning (i) equal interaction 
strength in all wires, q — qs, and (ii) isotropic tunnel- 
ing amplitudes, we have Qi — q and a — b. The RG 
equations read 



da 
dA ~ 


/?a = 


D- 


dc 
dA ~ 


/?c- 


~cl 


where 


A, 


— 



^{{q - a)(l - a)(l + 3a) + c^{2 + q + 3a)} , 

l)-3(a2 + c2)}. 

(13) 
By putting the /3- 



r/{l + 2g + 4a(q 
(a - qf + c^. 



functions equal to zero we find four fixed points, labeled 
iV, M, x"^ , C^ , at positions 



a-N = 


1, Cn = 







an = 


-1/3, CM 


= 


0, 


a^± = 


-1/2, c^± 


= 


±73/2 


n^+ — 


3-K^ 




"^+ - -1- 



(14) 



3(i^-l)2^ 



2Ky/K{K - 2) 
3(X-1)2 



The location of these points in the c — a-planc is shown in 
Fig. |3] We observe that all the FPs lie on the boundary 
curve of the domain of allowed conductances. We discard 
the unphysical solution a = q, c — 0, because |a| < 1 < 

\q\- 

Next we discuss the stability of these FPs. We start 
with the iV-point, located at a cusp of the boundary curve 
A, see Fig. 2. The cusp is infinitely sharp, leaving only 
a single direction of approach from the inside of A : a = 
1 — X , c ~ . Linearizing the RG equation for a in the 
small quantity x we find 
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FIG. 3: (Color online) RG flows for the fully symmetric chiral case, for different values of interaction. The four starting points 
are the same in all plots. The non-universal positions of the chiral C* points are shown by open circles. 



dx 
dK 



^2x{K-^ - 1). 



It follows that the A^-point is stable for repulsive inter- 
action, when K^^ — 1 > . 

We now turn to the x*-FPs, which are again located 
at a cusp of Curve A, Fig. [2j Therefore, again there is 
only one allowed direction for the trajectory towards or 
away from x '■ (oj c) = 5 (~-'^' ^ v3)(l — x) , where x is 
found to obey 



between the basins of attraction of the M, x+ and X- 
fixed points. So again, an 'economy' principle suggests 
the simple picture with only three stable fixed points", 
for 1 < i^ < 3. Our analysis shows, however, that this 
conjecture is not fully correct, as the new unstable fixed 
points C^ appear for K > 2 , allowing the point M to 
become stable. Examples of flow trajectories for interac- 
tion strengths K = 0.7; 1.5; 2.07; 2.7; 4; 15 are shown in 
Fig-i 



B. Chiral anisotropic case 



dx _ ^ ( '^K 

dK~ \iTlo 



-1 



We thus see that the x^-FPs are stable in the interaction 
regime 1 < K Ki . 

The M-point is located at a flat section of curve A and 
we therefore have a two-dimensional manifold of possible 
trajectories leading to it. Expanding the RG-equations 
in linear order around the Af-point we flnd 



da 
dK 
dc 
dA 



6; 



1- K 



K 



'^ i2 + Ky '■ 



It follows that the Af-point is stable for K > 2 . We note 
in passing that the M-point cannot be discussed within 
the Abelian bosonization approach, as the criterion a 



1 is not satisfied.— 



23 



M 



Whereas the TV, M, x^-FPs have been described be- 
fore, we find in addition a new pair of FPs C^ here. The 
location of these FPs is not geometrical, but depends 
on the interaction. We first observe that the C*-points 
are physically allowed (i.e. are located within the open 
domain A) only for sufficiently strong attractive interac- 
tion, K > 2 (they are always residing at the boundary 
curve A). At K — 2 they merge with the A/-point. As 
K increases the C* -points start to move away from the 
Af-point in opposite direction until they end at the x^~ 
points when K — 3 . Upon further increase of K the 
C* -points move beyond the x^ -points along the bound- 
ary of A , until they end at a = — 3, € = ±3, in the limit 
K —^ 00 . The stability analysis tells that the C^ -points 
are never stable. 

As a result, we conclude, that for K > 3 the only sta- 
ble fixed point is the M point, and not the unphysical D 
point a = —1, c = suggested iiP^. In Sec. 6.3 of that 
paper the authors say: "While the x± fixed points are 
stable against a small change in the flux, we do not know 
whether the M fixed point has such a stability. The sim- 
plest assumption is that it does not. In this case the RG 
flows might go to the x± fixed points for any non-zero 
flux (/), starting from arbitrarily small V. Alternatively, 
it is possible that the M fixed point is stable against 
adding a small flux. In that case, there would have to be 
additional unstable flxed points deflning the boundaries 



The above RG equations, Eqs. (12), describe the flow 
of a, b, c towards the stable fixed points (FPs) inside the 
body B of allowed conductance values (see Fig. IT]), de- 
pending on the initial conditions. The fixed points are 
defined by the zeros of the /3-functions on the r.h.s. of 
Eqs. (12). Some of the fixed points are outside the al- 
lowed region B in a — b — c-space and will be discarded as 
unphysical. There are altogether six physical FPs (see ta- 
ble ^ , which may be classified into two groups, universal 
(independent of the interaction) and non-universal FPs. 

The three universal FPs are labeled N,A,x^- These 
FPs may be thought to be of geometrical nature, as the 
corresponding S'-matrix elements are 1 or . The FP N 
describes the total separation of the wires, meaning zero 
conductance in all components. By contrast, A describes 
the separation of wire 3 from the perfectly conducting 
wires 1-2. Finally, x^ are the FPs generated by chirality 
(left- or right-handed). 

The non-universal FPs are labeled M,Q,C^. They 
only exist or are physical in a limited region of interac- 
tion constant space, mainly for attractive interaction. As 
shown in table |l] the conductances at these FPs depend 
in a complicated way on the interaction constants. 

The locations of the M and Q points have been given 
irP2l. We reproduce the results here in the present no- 
tation. Both points can simultaneously exist only for 
attractive interaction, 5,33 < 0. Notice that for equal 
interactions g — g^ we have q — qa = Qi, t = |1 -I- g| (see 
below) and we return to ( 14 ) . 

The FPs M and Q merge along a line in interaction 
constant space, defined by r = 0, where 



= yr 



g2 + 2Qi 



2, K. 



(15) 
4/3, and 



The end point of this line is at K 
was named "tricritical" in Ref. [22) . 

The location of the FPs C^ is best represented in terms 
of the angular variables 9,tp,S, as 



COS0 = --(2 + Qi) , cos?A = l, 

cosC = ±F^[3 - (2<Z - Qi)(2 + Qi)]"^ . 
5 + ^1 



(16) 



The C* points lie on the surface of the body of physically 
allowed conductances, since cosfA = 1 corresponds to the 



outmost points of our parametrization. The FPs C^ are 
only in the physical domain if the interaction constants 
satisfy the requirement | cos9\ < 1 , implying — 5 < Qi < 
1 , and < cos^ ^ < 1 leading to the bounds 

< 3 - (2g - gi)(Oi + 2) < 1(5 + Q^f 

The interval shrinks to zero when Qi — > —5, and the 
above inequalities demand q — >■ —3, which in turn gives 
93 — >■ —7. Remarkably, these values again correspond 
to K = 2 and K^ = 4/3, dubbed the "tricritical" point 
above. 

From the Table |T] one may also verify that the points 
C merge with x^ at 



Qi 



^ 



K- 



2A7I 



1 



(17) 



implying, in the symmetric situation, the condition K = 
^3 = 3. 



STABILITY OF FIXED POINTS 



The above set of equations ( 12 ) allows to perform a 



rather straightforward analysis of the stability of FPs. 
Assuming that we have a fixed point (ao, 60, cq), we may 
expand the RG equations in terms of the vector of small 
deviations x = {a — aQ,h — b^^c — Cq) = A(ai, &i, ci), to 
linear order in A ^ 1: 



d 
dK 



M 



x = 



(18) 



The stability of the fixed point is defined by the direction 
of the RG flow, which is determined by the eigenvalues 
Hi {i = 1,2,3) of the matrix M. All Hi are positive for 
the fully unstable FP, all jii are negative for stable FPs 
and in the the saddle point case (some Hi negative, some 
positive) the FP is again unstable. 

Two remarks are in order here. First, the matrix M is 
in general not symmetric, so that its eigenvectors are not 
orthogonal. Second, in the case of asymmetric interaction 
(9 7^ Qi) the expressions for M are not simple for the 
non-universal FPs, and should be analyzed numerically. 
Simplifications occur &t q — q^i — Qi, a& was already 
demonstrated for the M point in the non-chiral case irP^l. 

We mention that not all directions of the vector 
(ai,6i,ci) are permitted. This is because all FPs he 
on the surface of the body of physically allowed conduc- 
tances. Certain displacements x would take the point 
(a, 5, c) outside of the body. If the FP lies on a smooth 
part of the surface of the body, the requirement for x to 
be inside the body is not very restricting, because it only 
selects half of all possible directions. But in our case the 
situation is more complicated by the fact, that some FPs 
lie on ridges of the body. This situation happens with 
the universal FPs and is discussed below. 

In the following we present a detailed stability analysis 
for the different fixed points. The results are collected in 
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FIG. 4: (Color online) The RG phase portrait, showing re- 
gions with different set and character of the FPs. The region 
of attractive interactions is shown in more detail in the sec- 
ond panel. The numbers in the plot refer to its description 
in Table III] The three dots correspond to the special points 
K = K'i = 2, K ^ Ki = i, &nd K = 2, K'i ^ 4/3, discussed 
in the text. 



Table [TT| and in Fig. El showing the phase boundaries (i.e. 
the boundaries, separating different domains of stability) 
in the g — 53-plane. 

A corresponding figure for the non-chiral case (c — 
0) has been presented as Fig. 4 in our previous 
publicatiorP2l. By comparing the figures we observe that 
phases I- VI and 1,3 ',5 are characterized by FPs with 
c = and were already discussed in the non-chiral case. 
The new phases are labeled 2-4, 6, involving either the 
stable chiral FPs x* or the new unstable FPs, C*. 

It should be emphasized, that the scaling exponents 

obtained by our method for the "universal" FPs iV, A, 

X^ are in exact agreement with the results, found by 

I0] 



boundary conformal field theory^. Our results for the 
non-chiral FPs iV, A were reported earlier^. The par- 
ticular expressions for the x^ points below, Eq. (24), co- 



TABLE I: The positions of the fixed points of the RG equations, ( 12 1. The quantity r is defined in Eq. ( 15 1 



FP 
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b 
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TV 
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M 
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-1 
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|(kl -^)sign(g) 

|(|g|+r)sign(<?) 

\ (2g(0i + 2) - Q? + 1) 
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_i 

|((|g|-r)2-3) 

|((kl + ^)'-3) 

|(Qi(Qi+4) + l) 


t|(Qi- 


^l)\/3" 





+ 2 








c± 


-{2q- 


-Qi)(Qi + 2) 



TABLE II: Summary of the RG phase portrait, depicted in 
Fig. [4] The existence and stability of the RG fixed points is 
shown below. FPs A'", A, x always exist, and the stability is 
indicated by an "s" symbol. The stable and unstable FPs Af , 
Q, C axe denoted by "s" and "u", the "-" symbol is used if 
such a point does not exist. 
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incide with the exponents deduced from Table I of Ref. 
[20] for the scaling dimensions of the leading irrelevant 
boundary operators, A'*^. This becomes clear after the 
identification of our a^ with their 2 — 2A'<^, and after 
putting their gi = 52; 53 equal to our K; K^. 



A. N point 

In this case the matrix M is given by 

4(1 -Qi) Qi-q 

Mat = /at 4 - 3q - Qi 

3(2-g-Qi) 

In = [{q-l)(Qi-l)]-^. 



(19) 



Formally, there are three (left) eigenvectors (1,-1/3,0), 
(0, 1, 0), (0, 0, 1), corresponding to eigenvalues 4/(1 — g), 
l/(l-g)+3/(l-Qi), 3/(l-9) + 3/(l-Qi), respectively. 



However, as can be seen in Fig. [21 any displacement 
with ci 7^ is not allowed, as it would end outside the al- 
lowed domain. Therefore, when discussing the RG flows 
starting inside "the body", we should discard any dis- 
placements along c. We thus end with two exponents 
governing the RG flow towards N 



aN.i =4/(1 

1 

aN.2 



q) = -2{K-^-l) 



3 



1 



-{K- 



k; 



2) 



(20) 



q l-Ql 

which may be interpreted as the weak link exponents in 
the main wire and between the main wire and the tip, 
respectively (eq. (46) in^--). These values are negative 
for K, K3 < 1, so that the N point is stable in the case 
of repulsive interaction in all wires. 



B. A point 

The matrix M for the asymmetric A point is given by 



M, 



Ia 



4(Q 



1-1) Ql+q 







Qi-4-3q 







Oi 


-A' 


-3g 



fA = [iq + l)iQl 

We have three (left) 



1)] 



(21) 



eigenvectors (1,1/3,0), (0,1,0), 
(0,0, 1), corresponding to eigenvalues 4/(1 4- g), 1/(1 + 
q) + 3/(1 - Qi), 1/(1 + q) + 3/(1 - Qi), respectively. 
The last eigenvalue is doubly degenerate, and therefore 
we again have only two exponents 



aA,i = A/{l + q) = 2{l-K) 
1 3 

"A, 2 



1 



q 



1 -I 



(22) 



which may be interpreted as i) the weak impurity ex- 
ponent in the main wire and ii) the tunneling exponent 
in the main wire and the boundary exponent in the tip, 
respectively (eq. (47) in^^). At the "tricritical" point, 
K = 2, Kj, = 4/3, the second exponent vanishes. This 
is a manifestation of the fact, that the tricritical point 
is the endpoint of the line of stable fixed points in the 
interaction parameter space, see Fig. |4] 
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C. X point 

We consider only one of these chiral points, 
which we obtain 



for 



of the corresponding set of RG equations around these 
points leads to very cumbersome expressions for the ma- 
trix M, which we do not show here. Instead we list the 
three scaling exponents for the conductances. 
For M fixed point we have (for g < 0): 



4Qi +3q+f 
Mv = /x 3(gi + i) 
V3(2Qi - i) 
/x - 2[(2g+l)(2gi + l) 



2Qi+q 
Ql+ 6q 

V3{2q ' 

+ 3]-i. 



3 
2 

19 
2 



-V3(2gH 
-V3(20i 
9 



1) 
-1) 



(23) 



The three left eigenvectors are (-3, ^^^:^^, 2^3), 

(iSl^,^^^!^,-^), (^,V3,-2), corresponding 
to the eigenvalues 2/x(3g-Qi+4),4/x(Qi+2),3/p^(Qi + 
(7 + 4), respectively. The last eigenvector points out of 
"the body" and should be discarded, similarly to the sit- 
uation at the N point. 

The other two physical exponents are 



^x.i 



^X,2 



4(3g-Qi+4) 
(2g + l)(2Qi + l)+3 
o 8KK-S 

2K + Ky. + K^Kz 

8(Qi + 2) 



(24) 



(2g 
= 2- 



l)(2Qi + l)+3' 
AK{K + K:i) 
2K + Ka + K'^Ks ' 

For K3 = K we have both exponents equal 



a 



x,i 



*X:2 



2-8/^/(3- 



When the C^ and x* merge, which happens at Qi 
-2, the second exponent a^^2 in (24 1 vanishes. 



For K^ ^ K the form of the eigenvectors is not trans- 
parent, and it might be useful to restore the parametriza- 
tion in terms of Euler angles. The coordinates of the x^ 
points are then given hy 9 = vr/2, ^ = 0,7r, for arbitrary 
^. For x~, we put ?/; = and expand 6 = n/2 + 9', 
^ = TT + ^' , with 9',^' <C 1. Then, linearizing the result- 
ing RG equations, one finds a decoupled set of differential 
equations for 9',^' 



aj\/,i 



aM,2 



aj\/,3 = 



3r((g + r)^-9) 
2{q + T)i2q-T)^ 
3((g + r)^ + 3) 
' {q + T){2q~T) ' 
3(g + T + 3)(g2 



(26) 



2g - t2 + 2r - 3) 



2(g + T)(2g-T)2 



with the 

( q{{q+rf-9) 

defined in (15). 



corresponding 
-1,0), ((7 + T,-l,0), 



eigenvectors 
(0,0,1) and r 



For Q point we obtain (again, for g < 0) 
3r {{q - r)2 - 9) 



"Q,i 



"Q,2 



"Q,3 



2(g-T)(2q + r)2' 

3((g-r)^ + 3) 
{q-T)i2q + T) ' 
3{q-T + 3){q^ + 2q- 



(27) 



t2 - 2t - 3) 



2{q- 



T){2q- 
9((g- 



and the eigenvectors are ( g^,^ , iT^ w 
T,— 1,0), (0,0,1), correspondingly. 
For the C* point we get 



9) 



-1,0), {q 



12 



12 



ac,i 



"C,2(3) 



-3? + Qi + 
3 (3g - Qi 



2 gi 

f 4) (Qi 



6. 



5) 



2(Qi-l)(-39 

Qi 
Qi 



1 1 
-5^3 



39 + Qi + 2 


) 




V" 


/l , 8(Qi 


+ 2) (3(7 ■ 


-gi- 


-2) 


V ' (Qi 


- 1) (3? - 


-Qi + 


4) J 



(28) 



while the expressions for eigenvectors are too complicated 
to be presented here. One can check that at Qi — —2 the 
exponent ac,2 vanishes, it corr espo nds to C merging 



with X J a-s discussed after Eq. ( 24 ) 



d9' ^ , 2(3g-Qi+4) 

dK (2g+l)(2Qi + l)-^3' 

dj' ^ , 4(Qi + 2) 

dK (2q + l)(2Qi + l) + 3' 



(25) 



in correspondence with the above exponents for the con- 
ductance components. 



D. M, Q and C* points 

As seen in Table |lj the position of these FPs is not uni- 
versal, i.e. depends on the interaction. The linearization 



E. Lines of fixed points 

In our previous worWl^l on transport through a chiral 
y-junction at weak coupling we found that at special 
lines in the g — 53-plane lines of fixed points may appear. 
We conjectured that this may happen at any boundary 
in the g— (73-plane separating phases with different stable 
FPs, which still exist upon approach to the boundary. A 
precondition is that the two different stable fixed points 
should exist as separate fixed points at the respective 
boundary in the g — (73-plane, i.e. the FPs should not 
have merged or otherwise disappeared upon approach to 
the boundary. By inspecting the many phase boundaries 
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in Fig. 4 in that respect we find that only at two phase 
boundaries, (1) separating phases 3 and VI, at K = 1; 
K3 > 1 (i.e. g = 0, 33 < 0) and (2) separating phases 
1 and 3, at i(X + K'^) + K^^ -2 = 0, l<i^<2a 
Une of fixed points emerges. Examples of trajectories in 
that case, demonstrating the existence of a line of fixed 
points, have been given in our previous work for the weak 
coupling regime 
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VI. FULLY ANISOTROPIC Y-JUNCTION 

In the previous sections we presented a rather detailed 
analysis of the 1 — 2 - symmetric Y-junction in the pres- 
ence of magnetic flux. We now consider the effects of the 
absence of symmetry. As in previous cases, one might ex- 
pect that some of the stable fixed points found above will 
become unstable if the additional freedom of the RG-flow 
in the absence of the symmetry restriction is admitted. 
One may distinguish two ways in which the symmetry 
with respect to wires 1-2 may be broken. The first is 
provided by different properties of the wires 1,2, in par- 
ticular different interaction strengths 51 7^ 52- A sec- 
ond symmetry characterizes the junction, i.e. the way 
in which the third wire is attached. For instance, if the 



third wire is non-perpendicular to the main wire, one 
may expect some asymmetry in hopping to wire 1 and to 
wire 2. The resulting asymmetry in the S- matrix and the 
conductance matrix may grow in the process of renormal- 
ization. Such a scenario is implied in the work^". In order 
to compare our results with those irpH, we now discuss a 
possible anisotropy in the conductance tensor. 

In the following we confine ourselves for simplicity to 
the asymmetry of the second type, keeping in mind that 
the completely general case is fully covered by the gen- 
eral form of the RG-equations. The relevant formulas are 
listed in the Appendices KjlDl The S-matrix is now char- 
acterized by four Euler angles, and instead of Eq. (|6| we 
have four conductance parameters 



(29) 



with c 7^ — c. In the asymmetric non-chiral situation 
we have c = c, see^. Using the general KG equation 




(11) and formulas in Appendix PdI we obtain after some 



calculation: 



da 

dA 
db 

dA 
dc 
dA 
dc 
dA 



^ DY\a^i3b - 1 - 2Qi) + a(-3cc + g(l - b)) + (Qi - b){l + b) + c^~ ccq + <f} , 
= DY^{q{l - b){l + 26) + a{b - 1)(1 + 36 - Qi) ~ cc(2 + 36 + Qi)} , 
= Df^c (a(36 - 2Qi + 2) - q{2b + 1)) + c (-26 ^3c^ + Qi + l)}, 
= DY^{c (a(36 - 2Qi + 2) - q{2b + 1)) + c (-26 - 35^ + Qi + l)} . 



(30) 



with Di — (a — g)(6 — Qi) — cc. In the limit c = — c, we 
return back to ( [l2| . 

The investigation of the whole phase portrait defined 
by the set of equations, (301, is beyond the scope of 



the present study. For our purposes it suffices to check 
the stability of the already discussed FPs in the slightly 
asymmetric situation, which can be done as follows. Our 
previous analysis was confined to the surface c + c — 
in the space of conductance components a, 6, c, c. We 
now allow for a small asymmetry of the Y-junction, 
|c -t- c| <C 1 and expand the RG equations (301 to first 
order in h = c + c. As a result, we recover the three 



previously found equations ( 12 ) and the fourth equation 
takes the form: 



dh 
dA 



hD-\l-q + Qi+2,c^ 



-26(g + l)-a(2Qi-36-2)) 
(31) 



A negative (positive) coefficient of h on the right-hand 
side of the last equation corresponds to stability (insta- 
bility) of the previously determined FPs with respect to 
the asymmetric distortion of the S-matrix. 

It is clear, that if the FP does not exist in some region 
of interaction parameter space, or is unstable, then the 
additional check of the stability provided by (31) is not 



required. Only otherwise stable FPs should be subjected 
to this additional test. 

Performing this analysis, we find that the FPs collected 
in Table [n] are stable against asymmetry, with only one 
exception. Namely, the FP M is unstable with respect 
to the asymmetry perturbation /i 7^ in region VI. In 
all other situations the character of the listed FPs is un- 
changed. Particularly, the exponents for h around the 
points N and x are given by aN,2 and a^^2i Eq. (20) 



and (24), respectively. The RG flow from the M point 



in the region VI should apparently lead to one of two 
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FIG. 5: (Color online) The RG phase portrait, showing re- 
gions with different character of stable FPs. The numbers in 
the plot refer to the description in Table III] The three dots 
mark the points K — K3 — 2, K = K3 = 3, and K = 2, 
K2. = 4/3, discussed in the text. The three straight red lines 
correspond to {2K + K3)/3 — 2, 3, 4, respectively. The dashed 
line indicates the condition K = Kg. 



asymmetric FPs, different from A and discussed irP^l. 
These points A1.2 are defined by a = 1/2, b = —1/2, 
c = c = ±\/3/2 and correspond to the cases of either 
the first or the second wire detached from the remain- 
ing two wires. We will denote these points by A below. 
In addition to these, we also may have the counterparts 
of the point Q discussed above, at strong attraction. 
The unstable Q point separates the A and M points, 
as discussed irP^I. Two new asymmetric Qi,2 points ap- 
pear for strong repulsion and their position is given by 



a = i(2g(Qi 



2) 



n-i) 



b = i(Qi(gi+4) + l), 



c = c = ±1 (C^i - 1) ^{2q + Qi) (Qi + 2) + 3. 

We notice also that the asymmetric FP A cannot be 
affected by a finite value of h without violating unitarity, 
meaning that the scaling behavior is determined by only 
two exponents, Eq. (22). This situation is completely 
similar to what was discussed above, where the chiral 
perturbation was found impossible around the N point, 
and one of the perturbations impossible around the chiral 
X points. The corresponding boundary (deltoid) curve 
in the asymmetric non-chiral case, c = c, was shown in 
Fig. 1 of Ref.l^. 

Below we compare our results on the stability of differ- 
ent FPs with the results summarized in Fig. 5 of Ref.^ffl. 
In order to facilitate the comparison, we re-plot our Fig. 
[4] in terms of the Luttinger parameters, K, K3 and show 
the phase portrait in Fig. [5j 

We observe full correspondence of our results with 
those shown in panels (a), (b), (c), (d) in Fig. 5 of Ref.l^S 



for relatively small interaction strength. However, we 
only partly agree with results obtained in^'' for strong 
attraction, panels (e), (f), as discussed below. 

The three straight red lines in Fig. [5] designating the 
manifolds {2K + K^)/?, = 2, 3,4, correspond to the ver- 
tical lines drawn through the top vertex of the triangles 
shown in panels (d), (e), (f) in Fig. 5 of Ref.l^Sl, respec- 
tively. 

We see that the first line, {2K + i^3)/3 = 2, consecu- 
tively intersects the regions IV, V, VI, 3, 1, I, II, when 
one increases K from K — Q to K = 2. This sequence 
corresponds to the appearance of stable FPs in the fol- 
lowing order: N^ A, x^ ^ ^: in fuU agreement with Fig. 
5(d) of Ref:?2l. 

The situation is more involved at stronger attraction. 
The second line, {2K + ^3)/3 = 3, consecutively inter- 
sects the regions IV, V, VI, 3, 4, 3', 5, 1, I, II, in Fig. [5] 
upon the increase of K from K = Q to K = 4.5. From 
our Table [n] we see that this sequence corresponds to the 
FPs: N, A, x^. [X^+M), M,{M + A), A. It should be 
compared to the sequence TV, A, x^, A, A, shown in Fig. 
5(e) in^". We notice two differences here: one is related 
to the fact that Hou et al^^ do not include the M point 
in their analysis, therefore they do not show the combi- 
nation of several stable FPs, e.g. (x^+-A/) in region 4, 
cf. Fig. [3] at ii' = 2.07 above. Another difference con- 
cerns the region 3' which appears, e.g. si K — 2, K^ — 5. 
Our analysis shows that, given the initial symmetry with 
respect to wires 1 and 2, the RG does not flow away 
from this symmetry, and the stable FP in this case is M . 
This is in contradiction to the conclusion of the authors 
oPSl, who observe that the asymmetric points A are sta- 
ble and hence any RG flow should end there. In fact, one 



can show from Eqs. (30) that the symmetric point M is 

from the asymmetric points 
mentioned above. 



separated in the region 3' 
^1,2 by the unstable FPs (5i,2, 

The biggest difference occurs at even stronger attrac- 
tion. The third line, {2K + Kz)/i = 4, in our Fig. [5] in- 
tersects the regions IV, V, VI, 3, 4, 3', 6, 2, 1, I, II, these 
correspond again to the sequence of stable FPs : N , A, 
X^^ ix^+M), M, {M + A), A. This sequence should be 
compared to Fig. 5(f) irPSl^ which provides the list: TV, 
A, x^, A, D, A. We see that the bosonization analy- 
sis, apart from systematically ignoring the existence of 
the M point, proposes the appearance of the new fixed 
point D. The "Dirichlet" FP D violates the unitarity 
of the S-matrix, i.e. th e cons ervation of the number of 
particles; the authors o f^°*^° l explain that this property 
may be attributed to the formation of superconducting 
pairs close to the Y-junction. In particular, the D point 
is proposed as a stable FP even in the fully symmetric 
situation, K = K^. Our analysis does not indicate the 
existence of a FP with the properties of the D point, 
which lies outside the physically accessible region A in- 
dicated in Fig. [2] at the point a = 6=— l,c = 0. The 
relative position of the D point with respect to the body 
B' of allowed conductance values, available in the fluxless 
asymmetric case, is shown in Fig. l6] 



13 




context of superconducting fluctuation contributions. 
These are the Maki-Thompson and Aslamazov-Larkin 
correct iona^S. We showed, however, in our previous 



FIG. 6: (Color online) The position of the Dirichlet D point 
is shown in relation to the body of allowed conductances in 
the asymmetric non-chiral case. The fixed points N, Ai, M 
are also shown. 



It should be stressed again, that the scaling expo- 
nents at the points iV, x> A obtained irpSl coincide with 
those, calculated within our approach earlier and in the 
present paper. Let us look closer at the argument of the 
author&^° in favor of the D point. They observe that the 
pair hopping operator at the FP A becomes formally rel- 
evant for large attraction, which means that the RG flow 
generated by such a perturbation should move the sys- 
tem away from the A point. At the same time, they show 
that the perturbations around the proposed D points are 
irrelevant and the D point is thus stable. They do not ac- 
tually suggest any RG flow, connecting the FPs A and D, 
because these points are proposed by an ansatz, 1221 rather 
than determined as zeros of the beta-function of a set of 
RG equations. 

We observe further that the bosonization approach 
suggests the existence of the stable FP D at strong at- 
traction, whereas the existence of such a point is not 
excluded at weaker attraction, when it only becomes un- 
stable. Therefore it should be meaningful to discuss val- 
ues of conductance, which exceed the maximum allowed 
value for one-particle processes, for arbitrary interaction, 
in particular also for weak interaction. One may then 
hypothetically find examples of contributions from pair 
hopping, which may increase the conductance over its 
maximum value, even though the subsequent renormal- 
ization (i.e. higher order corrections) reduces this value 
to usual one. 

From a formal perspective, contributions to the Kubo 
formula of conductance involving fermion pair processes 
are certain vertex corrections, first discussed in the 



worlP^I that vertex corrections vanish at zero tempera- 
ture in the limit loL — > 0. It is known that the criterion 
ljL <C 1 is crucial for obtaining the correct value of the 
conductance of the clean Luttinger liquid with leada^nSl. 
This c riterion however does not appear in the analysis 
q j-io | 20 |^ where the presence of leads at a; > L is modeled at 
the end of the calculation by introducing "contact" resis- 
tances, which may be positive or negative depending on 
the interaction. After this modeling the authorJiSMl con- 
clude that the D point exists also in the presence of leads 
and is characterized by the above values a = b = —1, in 
violation of the unitarity of the single-particle S-matrix. 
The proposed picture raises a natural question about the 
role of the length of the interacting region L. In our 
fermionic analysis small L correspond to the absence of 
any renormalization. By contrast, the result in^*^ '^^ re- 
duces to the statement, that any finite region around 
the Y-junction should lead at strong attraction between 
fermions to the peculiar behavior that a fermion coming 
from one Fermi-liquid lead is scattered to another lead, 
and takes with it a partner fermion from a third lead. It 
is apparently suggested that the finite region close to the 
junction acts as a pump, adding extra charge from all 
leads to all leads. If valid, such a picture should have a 
precursor at moderate strength of interaction. However, 
we quoted above the absence of such contributions in the 
direct perturbative calculations up to the third order in 
the interaction. We conclude that the most plausible res- 
olution of the contradiction is that the ad hoc proc edure 
of modeling the non-interacting outer leads used by^^^^ 
is not correct. There is no reason why the "contact resis- 
tances" should be the same for the ideal (clean) T LL-wi re 
and wires connected by a junction, as assumed iip2Eol_ 



VII. CONCLUSION 

In this paper we presented results for the effect of 
fermion-fermion interaction on the transport properties 
of three quantum wires joined at a Y-junction. We 
employed the Luttinger liquid (LL) model of spinless 
fermions to calculate the S-matrix and the conductances 
of this system for the case of equal interaction constants g 
in the main wire and different interaction g^ in the third 
wire. Using a purely fermionic representation we are 
able to directly describe the physical situation of wires of 
length L adiabatically attached to reservoirs. We argue 
that in this case (at least at zero temperature) the scat- 
tering is purely elastic (excitation of particle-hole pairs is 
excluded) and therefore described by the single-particle 
S-matrix (the S-matrix is strongly renormalized by in- 
teraction effects). We calculate the renormalization of 
the S-matrix from a set of three renormalization group 
equations for three conductances. The RG-equations are 
derived by calculating the leading scale dependent con- 
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tributions to all orders in the interaction. These contri- 
butions are obtained by analytically solving a set of inte- 
gral equations (ladder summation). In our earlier work 
we have presented arguments suggesting that the selected 
contributions are the dominant ones, meaning that ad- 
ditional contributions do not change the scaling behav- 
ior near the fixed pointsi^. Additional support for the 
correctness of our method is coming from the fact that 
our results agree with known results in all cases where 
these results are well founded. The relative simplicity of 
our method allows, however, to consider more complex 
problems not addressed before, e.g. in the present case 
anisotropic Y-junctions at strong coupling. As demon- 
strated in our earlier work, the addition of even small 
non-symmetric perturbances to a highly symmetric sys- 
tem may change the RG-flow pattern in a qualitative way. 
For example, a second order contribution to the tunneling 
from a tip into a quantum wire is found to destabilize the 
fixed point corresponding to the tip separated from the 
ideally conducting wire^H. Similarly, a small asymmetry 
in the interaction constants of tip and wire of a chiral Y- 
junction leads to the appearance of lines of fixed points^*^. 
The fully symmetric chiral Y-junction has been studied in 
great detail by^Sl. These authors found the emergence of 
a pair of chiral FPs, x^ , which we confirm in our present 
analysis, extending the considerations to anisotropic in- 
teractions g ^ gs- They also conjectured a new and quite 
unusual fixed point termed D, which in their interpreta- 
tion involved "Andreev scattering", i.e. the tunneling of 
fermion pairs. Such processes necessarily involve particle- 
hole excitations, which we exclude from the beginning 
(these processes have vanishing phase space in the limit 
of zero excitation energy). The existence of the new FP 
was argued iniSl to be necessary to complete the RG flow 
pattern. We agree that new fixed points are necessary 
at large attractive interaction. However, we find in the 
isotropic case a set of additional unstable chiral FPs, C^, 
in the domain of attractive interaction, K > 2, serving 
a similar purpose. The emergence of C^ allows to stabi- 
lize the FP M. In the anisotropic case a further unstable 
FP, Q, is appearing. In a follow up tcJ^^I the properties 
of a fully asymmetric Y-junction have been considered 
recently^ni. Our detailed comparison with the results ob- 
tained by the latter authors reveals agreement in many 
aspects of this multifaceted problem. At strong attrac- 
tion, however, our results differ from theirs in two re- 
spects: we do not find their fixed point D, which violates 
the unitarity of the single particle S-matrix, but instead 
find additional fixed points not discovered by them. Fi- 
nally we like to point out that our formulation is rather 
general, allowing to derive RG-equations for junctions 
connecting any number of leads in any asymmetric way 
(see Eq.([TT])). The challenge in that case is to find a 
suitable parametrization of the S-matrix. Work in this 
direction for the four-lead junction is in progress. 
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Appendix A: S-matrix parametrization and 
re-phasing 

The traceless Gell-Mann matrices, Xj, with j = 1, ... 8 
are the generators of the SU(3) group and are listed 
elsewhere.^"* Together with the unit matrix Aq, they form 
a basis for representing the 5-matrix of a three-lead 
junction. We note the property Tr[AjAfc] = 2(5jfe, with 
J, fc = 0, . . . 8. 

Up to an overall phase, an arbitrary matrix S, belong- 
ing to representations of the SU(3) group, may be defined 
by its Euler angles parametrization as foUows^S : 



g — Ugi>^aSjj^iX8a4,(V3/2) 



u 



JX3a2/2 i\2i/2i\3ip/2 



(Al) 



JJ ^ giA30/2gjA2«/2giA3Q3/2^ 



While this expression for S contains eight Euler angles, 
for our purposes it is convenient to introduce an addi- 
tional degree of freedom by adding a ninth angle and 
defining 



5";^ — g^-^8"l(v^/2)5'_ 



(A2) 



One may now use the redundancy in the latter repre- 
sentation to derive a more convenient representation of 
S. Using symbolic computer calculations, one can verify 
that Si is invariant under the following change 



V" ^ V' + 3/3, 0^^-3/3. 



(A3) 



for arbitrary f3. This means that we may set, e.g., (jj — 0, 
by choosing /? = 0/3. Doing this, we arrive at an equiva- 
lent representation of the S'-matrix in terms of eight pa- 
rameters. This representation contains four re-phasing 
angles ai, . . . a^. These do not appear in the expressions 
for the conductances. In the main text we therefore dis- 
carded the re-phasing angles, which leads to a simpler 
parametrization of the S'-matrix. In sections Il]|IV[ we 
set ^ = TT — ^ and work with the three angles 9, ip, ^. 
In sections |VJ |VI[ Appendix |D] we also refer to the most 
general case with four parameters, 0, ip, ^, ^. 
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Appendix B: Correspondence with other authors 

Let us first establish the correspondence with the no- 
tation of Oshikawa et all^ in their Appendix A. They 
obtain the S'-matrix in the form 



S = e~'^''' [1 - e^'TM] ^ [1 - e'^'TM] , 



(Bl) 



where k is the wave- vector, and F is the amphtude in the 
tight-binding Hamiltonian, describing the hopping from 
one wire to another. The matrix M is given by M = 

e^^/s^t _,_ e-"^/^B, and B 



0, 1, 

0, 0, 1 

1, 0, 



Here 6 measures 



chirahty, i.e. is a normalized flux through the Y-junction. 

For small F we have (up to an overall phase) S ~ 
1 + 2i sin k TAI. In the tight-binding, particle-hole sym- 
metric case, considered below, we have k ^ kp = 7r/2 
and 5 ~ 1 -h 2iTM. 

We can express M in terms of Gell-Mann matrices A^ 
as 

M = cos |(Ai + A4 + As) + sin |(A2 - A5 + A7) . (B2) 

It is natural to parametrize the ^-matrix in a mini- 
mal way by the expression S = exp [z^M] , where 9 — 
2F in the limit F ^ 1. It turns out, however, that 
the expressions for conductances are simplified using a 
slightly different parametrization. Instead of the pair 
(6'cos(0/3), 6'sin(0/3)) in Eq. ([B2]) it is better to use the 
pair (01, 62) with 



6'cos(0/3) =6'i/3, 



'sm 



(0/3)- 02/^3. 



After this convention, we arrive at the S-matrix in the 
form 

S = exp (|0i(Ai + A4 + Ag) + ;^^2(A2 - As + Ay)) 
= e-'^^'^{r + t+B + t^B^) 

r, t+, t. 

t+ I (B3) 



e-^'^'^ \t , r, U 



t, t_ 



with 



r= i(e*^i+2cos02), 
t± = \{e'^' - cos6'2 ± \/3sin6'2) . 



(B4) 



The connection with parametrization (Bl I is established 
by taking Oi,2 <C 1 : 



6*1 = 3Fcos(<?!)/3), 62 = y3Fsin(0/3), 



(B5) 



here F denotes the hopping amplitude, and = 7rn is 
standing for the flux, piercing the junction, with n the 
number of flux quanta. We have in special cases 



3F, 02-0, 

|^2|-3F/2, 



mod 3 
: ±1 mod 3 



(B6) 



The components of conductance, a = b, c, Eq. (pi, are 
given by 



||r|2- i = i(cos26l2-l-2cos6liCos6'2) 



Va 



i\U 



t\ 



sin 02 (cos 01 — cos 02) 



(B7) 



We note that a and c reach their extremal values at 
cos 01 = ±1, and these values describe the boundary of a 
domain of the form 



a 



= i(cos202 + 2cos02), c= |sin02(l-cos02), (B8) 

This curve is called deltoid, and is depicted in Fig. [2] 
The Abelian bosonization approach may be applied to 
the problem when the condition a^ -I- c^ = 1 is satisfiecP^l. 
One observes that there are only three such points in Fig. 
[2] the corners of the deltoid. 

The simple relation (B5 1 holds only for small F. For fi- 



nite values of F we calculate the conductance components 



directly from (Bl| and obtain 



12F2(1 + F2) 



(1 + 3F2)2 + 4F6 cos2 



8v^F3 



sin( 



(B9) 



(l + 3F2)2+4F6cos2( 



These expressions correspond to Eq. (91) irP^ after a 
change F — > t. The quantities (a, c) lie within a deltoid 
curve for arbitrary F g (0, 00), (j) € (0, 27r), as shown in 
Fig. [2] The proper deltoid curve is given by Eqs. (B9) 
at = 7r/2 and F e (0, 00) 

It was argued in^" that the free-fermion description 
of the S'-matrix is not fully applicable in the interacting 
case. Instead, it was proposed to use the bosonization 
approach and to find an analog of small F in the bosonic 
Hamiltonian. When abandoning the fermionic formula- 
tion and working directly with currents (densities), one 
does not find restrictions on the values of a, c, except 
for the condition a- + c < 1, mentioned above. It was 
hence implied irP^that such restrictions, imposed by the 
unitarity of S'-matrix, are relaxed in the interacting case. 

For completeness, we quote here the lowest-order RG 
equations in terms of 0i 2 



d n 

d n 

5A^2 - 



'1 COS 02, 

ig sin 02 (cos 01 



-g sin 

1 



2 COS 02). 



(BIO) 



The resulting RG flows are shown qualitatively in the two 
first panels in Fig. l3J obtained for moderate interaction 
strength, K — 0.7, 1.5. 

The point of maximally open Y-junction, or M-point, 
is given by a = —1/3, c = and is obtained in two cases. 
First, in the absence of the flux, (p ~ 0, and F ~ 1/a/2, 
and second in the case of maximum flux (j> — ±7r/2, and 
infinitely strong hopping, F — >■ cxd. It is worth to com- 
pare this observation with Fig. 12 of Ref.'ISi^ where the 
case, 4> — ±7r/2, F — >■ 00, was associated not with the M 
point, but rather with a different FP labeled D. In our 
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FIG. 7: Feynman diagrams depicting the integral equations 
for the renormahzed interaction, Eqs.(C2l and (CI I. A sim- 



ple wavy hne stands for the initial interaction matrix g, the 
double wavy line defines the quantity L, leading to the effec- 
tive interaction strength in the RG equations for the conduc- 
tances. 



analysis we showed that the M point is the only stable 
FP for ii' > 3, as indicated in Fig. [3] above. Taken from 
this perspective, our results for K > i do not contradict 
the results obtained irfiSl^ since both = 0, F = l/-\/2 
and (/) — ±7r/2, F — >■ oo correspond to the same FP, the 
M point. The important distinction is, however, that 
our analysis describes the renormalization of the conduc- 
tances fully in terms of the one-particle S'-matrix, where 
the D point would be in the unphysical regime. By con- 
trast, we find the new FPs C*, not discussed irPSl. 



Appendix C: RG equation for conductances 

In this section we sketch the derivation of the formula 
([8|. This derivation closely follows our papeil^Sl, dealing 
with a junction of two non-equal Luttinger liquid wires. 
The formalism presented in that paper can be naturally 
generalized to an arbitrary number of wires, connected 
by a single junction. 

In the lowest order of perturbation theory the renor- 
malization of the S-matrix and the conductance matrix 
is proportional to the strength of interaction, Eq. (l7|. 
In higher orders of interaction the interaction constant, 
entering the RG equation for the conductance, is an effec- 
tive quantity, described by a function which is non-local 
and includes retardation effects. 

As shown irP^the leading terms in perturbation theory 
form a ladder series, which may be summed up analyti- 
cally to give a function L. The latter obeys the integral 
equation: 



L(a;,j/;a;) 
L2(a;,y;a;) 



2TTg,5[x - y) 




y gn(z-a:,a;), 

with the fermionic loop n(x,a;„) = {2'k)^^[5{x) — 
|w„|0(xc^„)e---"). 

It can be shown, that in the case with different strength 
of interaction in different semi-wires, it is more conve- 
nient to solve this equation in two steps. This is graphi- 
cally depicted in Fig. [7j whereas the details are presented 
elsewhere .1^ 

At the first step, the auxiliary quantity C is introduced 
as L, taken at Y = 0. It satisfies the equation 



C{x,y;uj) = 2TT^5{x-y) 



-'Uj{x-\-z) 



iwe 



-cj| j; — s| 



dz 



C{z,y;uj). 



(CI) 



1 



g?. The above condition 



with gj = Sj /d3 and dj - 

Y = is not intuitively evident, and it describes the 
absence of correlations between the densities of incoming 
and outgoing electrons. In this special case, the effects 
of scattering are absent and the effects of the interaction 
amount to the usual dressing of the incoming density by 
the outgoing one. Strictly speaking, the condition Y = 
cannot be realized for any S-matrix, as it violates the 
charge conservation, and should be viewed as vaguely 
corresponding to the situation "in the bulk" , far away 
from the junction. 



Solving Eq. (CI), we then obtain the renormalization 



of the plasmon velocity and certain prefactors of the in- 
teraction strength in the individual wires, which eventu- 
ally become Luttinger coefficients. 

At the second step, we use the expression for C thus 
obtained in the equation for L, which now includes the 
knowledge about the scattering Y : 



(C2) 



L{x,y;uj) ^C{x,y;uj) + -- dzidz2 

^TT Jl 

X C{x,zi;uj)Ye-^^'^+''^L{z2,y;uj), 



The solution of the above integral equation ( C2 ) is rela- 



tively simple. Substituting the result into the expression 
for the conductances one observes that retardation and 
non-locality effects, contained in L, disappear in the limit 
LuL — )• 0. We then arrive at the RG equation containing 
the expression for the effective interaction in the form 

g 

main text of the paper. 



ladder ^ 2(Q-Y)-i. This is Eq. M quoted in the 



Appendix D: Fully asymmetric Y-junction 



In the most general case the S matrix, (Al), is de- 
fined by eight parameters, but only four parameters, 9, 
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ip, ^, ^ define the conductance. The parameters of the 
conductance matrix ( 29 ) are given by 



cos"' e+1 



cos ^ cos £, — cos 9 cos tp sin ^ sin ^ , 



6= i(3cos2 6'-l), 



(Dl) 



-_V3 



COS £_ sin^ , c 



-_V3 



COS 1^ sin^ 6' . 



-Tr 



The general equation (111 contains a quantity with four 
indices, Fjkim 



w,1w/l 



This quantity is en- 



tirely defined by the parameters of ^'-matrix, as opposed 
to the interaction, entering the definition of g^; "^^^ . It 
is clear from the properties of Gell-Mann matrices, that 
Fjkim = if any of its indices equals 3; hence there are 
only 2^ = 16 non-zero components of Fjkim ■ An analy- 
sis similar to the one presented in the main text for the 
1 — 2 symmetric junction shows that Fj^im can be ex- 
pressed through the conductance components, a, b, c, c 
in a compact way. In view of definition ( 29 1 it is con- 



venient to refer to pairwise combinations of indices as 



{11} = a, {12} = c, {21} = c, {22} = b, so that i^nai 
is now denoted as Fas ', Feb stands for i^i222 etc. In this 
notation we have 



i^bfc = 1 + 6 - 26^ , 
F,, = l-b~2c\ Fsc = l-b- 2(? , 

Fab = Fba = Fee = Fee = a(l - 6) - CC , 

^ ac ^ ca C ZQC , ^ ac ^~ ^ ca ^ ZdC , 

Fbe = Feb^-i2b+l)c, 

Fbe = Feb^-{2b+l)c, 



(D2) 



These expressions stem from the structure of the SU(3) 
group, and the reduction to the SU(2) group is obtained 
by setting b — 1, c = c — 0. The latter choice leaves 
only one non-zero element Faa — 2 ~ 2a} in the above 
formulas and we return to the case of a junction between 
two Luttinger liquids considered earlier. 
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